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Abstract. In this contribution we summarize recent progress in understanding the shear viscosity of strongly correlated 
dilute Fermi gases. We discuss predictions from kinetic theory, and show how these predictions can be tested using recent 
experimental data on elliptic flow. We find agreement between theory and experiments in the high temperature regime T:$>Tp, 
where Tf is the the temperature where quantum degeneracy effects become important. In the low temperature regime, T ^Tp, 
the strongest constraints on the shear viscosity come from experimental studies of the damping of collective modes. These 
experiments indicate that rj/s <0.5ti/kB, where 7] is the shear viscosity and s is the entropy density. 
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INTRODUCTION 

Experiments carried out at the relativistic heavy ion col- 
lider (RHIC) have demonstrated that the quark gluon 
plasma is a very good fluid [1, 2, 3]. A quantitative anal- 
ysis of the observed elliptic flow in the framework of vis- 
cous relativistic hydrodynamics shows that [4, 5] 



^ < 0.5 ' 



(1) 



where tj is the shear viscosity and s is the entropy den- 
sity. The best fit to the data corresponds to even smaller 
values, r]/s :^ (0.1 —Q.2)h/kB. This number is close to 
a proposed lower bound, rj/s = h/ {Anks), on the ratio 
of shear viscosity to entropy density that is saturated in 
the strong coupling limit of a large class of field theories 
that can be analyzed using the AdS/CFT correspondence 
[6, 7]. 

The RHIC results raise the question whether nearly 
perfect fluidity is a phenomenon that is specific to rel- 
ativistic gauge theories, or whether it is a more general 
effect that also appears in other strongly correlated quan- 
tum fluids. If that is the case we may also ask whether 
there are universal aspects of nearly perfect fluidity that 
go beyond the relation Tj/s ~ /z/ (A-nks). Other universal 
features could include bounds on other transport coeffi- 
cients hke diffusion constants or relaxation times, or con- 
straints on spectral properties of the theory. In the context 
of interpreting the RHIC results we are particularly in- 
terested in the question whether nearly perfect fluidity is 
consistent with a quasi-particle description of the fluid. 

Cold, dilute Fermi gases in which the interaction be- 
tween the atoms can be tuned using an external magnetic 
field provide a new paradigm for strongly correlated 
quantum fluids [8, 9, 10]. These systems have been real- 
ized experimentally using optically trapped alkali atoms 
such as ^Li and '^'^K. These two atoms are fermions be- 



cause they posses is a single valence electron and the nu- 
clear spin is integer. When a dilute gas of alkali atoms is 
cooled to very low temperature, we can view the atoms 
as pointlike particles interacting via interatomic poten- 
tials which depend on the hyperfine quantum numbers 
of the valence electron. A Feshbach resonance arises if 
a molecular bound state in a "closed" hyperfine channel 
crosses the threshold of an energetically lower "open" 
channel. Because the magnetic moments of the states in 
the open and closed channel are in general different, Fes- 
hbach resonances can be tuned using an applied magnetic 
field. At resonance the two-body scattering length in the 
open channel diverges, and the cross section o is hmited 
only by unitarity, a{k) =47i/k^ where k is the relative 
momentum. 

In the unitarity limit details of the microscopic inter- 
action are irrelevant, and the system displays universal 
properties. The dilute Fermi gas can be described by the 
effective lagrangian 



^ = xi/Uido + 



2m 



V- 



(2) 



were \j/ is a two-component fermion field with mass 
m. The coupUng constant Q is related to the i-wave 
scattering length a. The precise form of this relation 
depends on the regularization scheme. In dimensional 
regularization one finds Co = {4na)/m. Terms with more 
derivatives or higher powers of the field are related to 
effective range corrections, higher partial waves, and 
many-body forces. All of these terms are irrelevant in the 
universal hmit. The two-body scattering cross section is 



aik) 



1+a^k^ ■ 



(3) 



In the Umit a «> the theory is parameter free, strongly 
interacting and scale as weU as conformally invariant. 



THERMODYNAMICS 

In the weak coupling limit a the equation of state 
of the dilute Fermi gas is that of a free two-component 
Fermi gas, 



Po{nJ) = -2TX^iLi„2{-t;-'), 



(4) 



where Xhb = [(27r)/(m7')]'/-^ is the de BrogUe wave 
length, Lia{x) is the Polylogarithm function, and C, = 
cxp{—jj./T) is the fugacity. At unitarity scale invariance 
implies that the equation of state is of the form 



(5) 



where h{^) is a universal function. This function can be 
calculated using the Virial expansion at high tempera- 
ture, C ^ 1, but it is a non-perturbative quantity at low 
temperature. The equation of state is known to about 
10% from experimental measurements [11, 12, 13] and 
quantum Monte Carlo simulations [14, 15, 16]. At zero 
temperature h{0)/2 = E, where £, ~ 0.4 is known as 
the Bertsch parameter. Scale invariance also impUes that 



3 ' 



(6) 



where e is the energy density. At low temperature the 
attractive interaction between the fermions leads to su- 
perfluidity. The critical temperature is Tc c^: O.lSEp [15], 
where Ef is the Fermi energy. The Fermi energy of the 
interacting gas is defined by Ep = kp/ {2m) where kp = 
(37r^n)^/^ is the Fermi momentum and n the fermion 
density. 



HYDRODYNAMICS 

At large distances and long times deviations from equi- 
librium are described by hydrodynamics. For simphc- 
ity we will consider the unitary Fermi gas in the nor- 
mal phase. In that case there are five hydrodynamic vari- 
ables, the mass density p = mn, the flow velocity v, and 
the energy density ^. These variables satisfy five hydro- 
dynamic equations, the continuity equation, the Navier- 
Stokes equation, and the equation of energy conserva- 
tion, 

dipVi) 
dt 

dS 



Vji^ = 



0, 


(7) 


0, 


(8) 


0. 


(9) 



The total energy density is the sum of the internal energy 
density and kinetic energy density, = e+ \pv^. These 



equations close once we supply constitutive relations for 
the stress tensor 11,, and the energy current as well as 
an equation of state. As explained in the previous section 
the equation of state is P = jE. The stress tensor is given 
by 

Ilij=pv,Vj+P8,j + 8llij, (10) 

where 5Ylij is the dissipative part. The dissipative contri- 
bution to the stress tensor is 511,^ = —^aij with 

..^(v,.. + v,.-?,xv...)), (.1) 

where 77 is the shear viscosity and we have used the fact 
that the bulk viscosity of the unitary Fermi gas is zero 
[17]. The energy current is 



.li 



(12) 



where w = e + /" is the enthalpy density. The dissipative 
energy current is 



5jt = 5nijVj-KViT, 



(13) 



where T is the temperature and K is the thermal con- 
ductivity. We note that the temperature T = T{n,P) is a 
function of the density n = p/m and the pressure. In or- 
der to determine T we need the equation of state in the 
form P = P(n,T). Universality implies that P{n,T) = 
m~^rr'l^ f{mT /n^l^) where /(x) is a universal function 
that is related to the function h{Q defined in equ. (5). 
The situation simplifies in the high temperature limit 
where P = nT. Universality also restrict the dependence 
of the shear viscosity and thermal conductivity on the 
density and the temperature. We can write 



K{n,T) 



mT 

(niT \ n 
J m 



(14) 
(15) 



where a„{y) and (Jniy) are universal functions of >' = 
mT /n^l^. The relative importance of thermal and mo- 
mentum diffusion can be characterized in terms of a di- 
mensionless ratio known an the Prandtl number, Pr = 
Cprj/{pK), where Cp is the specific heat at constant pres- 
sure. In the high temperature limit Cp = p/m and Pr = 
Un/On. Kinetic theory predicts that in this limit Pr = 2/3 
[18]. 



KINETIC THEORY 

Near Tc the transport coefficients tj and K are non- 
perturbative quantities that have to be extracted from ex- 
periment or computed in quantum Monte Carlo calcula- 
tions. At high temperature (and at very low temperature. 



T -^Tc, see [19]) transport coefficients can be computed 
in kinetic theory. The shear viscosity was first computed 
in [20]. Here we will follow the recent work [21] which 
also considers the frequency dependence of the shear vis- 
cosity. 

In kinetic theory the stress tensor is given by 



n,, 



d p p'p-' 
(2;r)3 m 



fp, 



(16) 



where fp — f{t,x,p) is the distribution function of 
fermion quasi-particles and the factor 2 is the spin de- 
generacy. The equation of motion for fp is the Boltzmann 
equation. In order to extract the shear viscosity it is use- 
ful to consider the Boltzmann equation in a background 
gravitational field. In this setting correlation functions of 
the stress tensor can be determined by computing varia- 
tional derivatives with respect to the background metric. 
The non-relativistic limit of the Boltzmann equation in a 
curved space chai^acterized by the metric gij is 



(- 



p' d 
m dx' 



) dp' 



/(r,x,p)=C[/], (17) 



where F^^ is the Christoffel symbol and C[f] is 
the collision integral. We consider small devia- 
tions from equilibrium and write f — fo + 5f with 
/o(p) = Mp'P^8ij/{2mT)). We also write gij = 5ij + hij 
and linearize in hjj and 5f. We get 



dt 



p' d 
m dx' 



/o(l-/o) 
2mT 



p'p'kj^C[5f]. (18) 



This equation can be solved by making an ansatz for 5f. 
We go to Fourier space and and write 



5/(co,^,p) = C(j/o(l-/o) 



p'p' 



2mT (0 



le 



(19) 



where B,^-^ — B, ^-^{cOjk), hjj^ are the traceless/trace parts 
of hij, and v'p= p'/mis the quasi-particle velocity. Insert- 
ing this ansatz into the Boltzmann equation we can solve 
for ^^'^ and then compute 5f and T/j. Matching the re- 
sult to hydrodynamics determines the shear viscosity. In 
the limit A: — > this can be done analytically. The zero 
frequency limit of the shear viscosity is 



15(mr)3/2 



1 



32V7r \ l/{3mTa^) a^O 



(20) 



We observe that the shear viscosity is large in the weak 
coupling limit a 0, and that 77 decreases as the scatter- 
ing length is increased. We note, however, that the shear 
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FIGURE 1. Time evolution of the aspect ratio As of a de- 
formed cloud as a function of the dimensionless time variable 
t = (Bx?, where (Bx is the harmonic oscillator constant of the 
transverse confining potential. The solid lines show the analyt- 
ical result for the ideal evolution, the dashed lines correspond 
to the dissipative solution without reheating, and the dotted line 
shows the solution of equ. (26-29). The points are from a nu- 
merical solution of the hydrodynamic equations. The viscosity 
parameter is j8 = 0.066, see equ. (30). 



viscosity saturates when the scattering length becomes 
comparable to the de Broglie wave length a ^ X^g ^ 
{mT)^^l'^. In this limit r\ only depends on Xds but not 
on the density or the scattering length. The frequency de- 
pendence of the shear viscosity is 



77(0)) = 



I? 



1 



(21) 



where T = (3T7)/(2e) is the relaxation time [22]. The 
relaxation time controls the time scale over which the 
viscous stress tensor relaxes to the Navier-Stokes form 
given in equ. (11). We observe that relaxation is fast in 
the limit where the viscosity is small. We also observe 
that the viscosity satisfies a sum rule which only depends 
on thermodynamic quantities. 



^ / d(or\{(o) = I 

TT Jo 3 



(22) 



A modified version of this sum rule which contains an 
extra short time (high frequency) contribution can be 
proven in the full quantum theory [23]. 



ELLIPTIC FLOW 

The first experiment that demonstrated nearly perfect flu- 
idity in the dilute Fermi gas was the observation of ellip- 
tic flow by O'Hara et al. [24]. The experiment involves 
releasing the Fermi gas from a deformed, cylindrically 



symmetric, trap. The density evolves as 



n(xx,x,,f) : 



1 



no(jc±^±(0)^A(0) 7 (23) 



where the transverse and longitudinal coordi- 

nate, b±{t),b^{t) are scale factors, and nQ{x±,x^) is the 
equilibrium density of the trapped system. The initial 
system is strongly deformed, A/;(0) = ^ 
1 . Hydrodynamic evolution converts the large transverse 
pressure gradient into transverse flow. As a consequence 
the aspect ratio AK(f) grows with time and eventually be- 
comes larger than one, see Figs. 1,2. 

Viscosity slows down the transverse expansion of the 
system. In order to quantify the effect of shear viscosity 
we have to solve the Navier-Stokes equation for the 
expanding cloud [25, 26]. The Navier-Stokes equation is 



dt 



+ V-V v,=/,+ 



(24) 



where fi — (V,f)/;7i is the force. With the help of the 
Navier-Stokes equation the energy equation can be writ- 
ten as 



d_ 
dt 



-vV- 



fi 



+ (V,-V;)/;-^(V,V,-V,)^ 



2V^ 

3 n ' 



(25) 



where q — 



Oij) is the heating rate. In general these 



equations have to be solved numerically on a finite grid. 
In [26] we showed that under certain assumptions a 
very accurate scaling solution can be derived. We will 
assume that the local shear viscosity is proportional to 
the density, r\ — a„n, where a„ is a constant. We will 
also assume that the systems remains isothermal and that 
heat conductivity is not important. 

The basic idea is to postulate that the velocity field and 
the force are linear in the coordinates. If the velocity is 
linear and rj ^ n then all terms in equ. (24) are linear in 
Xi. Also, equ. (25) is independent of the pressure and all 
the remaining terms are linear in Xi. We write /, = fl,x,, 
V,- = a,x, (no sum over /) and use the scaling ansatz 
(23) for the density. The continuity equation requires 
a,- = bi/bi. The scale parameters a, and bj are determined 
by the coupled equations 
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FIGURE 2. Data for the aspect ratio versus time, from [27]: 
Top Black, E = 0.6Ef- Pink, E = 23Ef; Blue, E = 33Ef ; 
Bottom Red, E = A.6Ef. Solid curves: Hydrodynamic theory 
with the viscosity as the fit parameter. 
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(29) 



where Ox , CO, are the oscillator frequencies of the har- 
monic confinement potential (before the gas is released). 
The parameter /3 is defined by 



i3 = 



a„ 



1 



(3A^l)i/3 (£o/£f)' 



(30) 



where is the number of atoms, A = A^(0) the initial 
aspect ratio, and Eq/Ef the initial energy in units of 
Ef = {3Niy/^NC0j_. The initial conditions are bj_{0) = 
b,{0) = 1, ^x(O) = k{0) = 0, and flx(O) = col, «z(0) = 
col 

Dissipative effects fall into two categories. The terms 
proportional to j3 in equ. (26,27) correspond to friction 
- shear viscosity slows down the expansion in the trans- 
verse direction. The dissipative terms in equ. (28,29) de- 
scribe reheating - shear viscosity converts some kinetic 
energy to heat which increases the pressure and even- 
tually re-accelerates the system. The interplay between 
these two effects can be seen in Fig. 1. As expected, fric- 
tion slows down the growth of Ak(0 as compared to the 
ideal evolution. Reheating reduces this effect by about a 
factor of 2. We also note that the solution of the equ. (26- 
29) agrees with very well with numerical solutions on a 
three dimensional grid. 

In order to make comparisons with data we have to 
take into account that a„ is not a constant. In [25, 26] 
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FIGURE 3. Trap average {as) = {ri/s) extracted from the 
damping of the radial breathing mode. The data points were 
obtained using equ. (33) to analyze the data published by Kinast 
et al. [28]. The thermodynamic quantities (S/N) and E{)/Ef 
were taken from [11]. The solid red and blue lines show the 
expected low and high temperature limits. Both theory curves 
include relaxation time effects. 



we argued that a„ in equ. (30) should be interpreted as 
the trap average of the local ratio of shear viscosity over 
density. 




In the dilute corona of the cloud the local shear viscosity 
is independent of the density and equ. (31) is not well 
defined. This problem can be addressed by taking into 
account that the viscous stresses relax to the Navier- 
Stokes value on a time scale T that becomes large as the 
density goes to zero, see equ. (21). A relaxation model 
for (a„) was studied in [25]. An even simpler model can 
be constructed based on the assumption that the shear 
viscosity relaxes to its equilibrium value at the center of 
the trap and is proportional to the density in the dilute 
corona. This implies that r\{x) ~ T7(0)(n(x)/n(0)). This 
parametrization agrees with the relaxation model at the 
30% level. It was used by Cao et al. [27] to analyze the 
data shown in Fig. 2. The hydrodynamic curves shown in 
Fig. 2 were obtained with Tj — riQ{mT)^^^ and Tjo 0.33. 
This agrees quite well with the prediction of kinetic 
theory Tjo = 15/(32 V^) ~ 0.26, see equ. (20). 

COLLECTIVE MODES 

At temperatures T < Tp dissipative effects on elliptic 
flow are very small. In this regime accurate measure- 
ments of the shear viscosity can be obtained by analyz- 
ing the damping of collective modes [28, 29, 30]. The 
analysis presented in the previous section can easily be 
extended to that case. A solution that describes a radial 



breathing mode is given by 

bj_{t) = 1 +flxcos((»f)exp(-j3(»xf) ■ (32) 

where a± <^ lis the amplitude, o = (10/3)'/^ Ox is the 
frequency, and j3 is the parameter defined in equ. (30). 
The experimentally measured damping rate F can be 
used to estimate (a„). We find [25] 

(a.) = (3«)'/3(i:)(|i). (33) 

In order to compare with the heavy ion data and the 
proposed string theory bound it is interesting to convert 
the ratio rj jn to the ratio r\ js of shear viscosity to entropy 
density. This can be done using measurements of entropy 
per particle published in [11]. In Fig. 3 we show an 
analysis of the collective mode data obtained by Kinast 
et al. [28] using equ. (33). The solid lines show the 
prediction of kinetic theory in the limits T <C 7f and 
T ^ 7>. We observe that 77/5 is consistent with kinetic 
theory for T >^Q.5Tf. We also find that in this regime 
the data from collective modes are in agreement with the 
results based on the elliptic flow data [27]. The data do 
not show the expected behavior at very low temperature. 
We should note, however, that for T <C 7> the mean free 
path is very large and dissipative hydrodynamics is not 
applicable. 

We observe that rj/s reaches a minimum close to 
the phase transition temperature, and that the minimum 
value is Tj/i ~ 0.5 (in units of Ti — kg — 1). We emphasize 
that this result refers to trap averaged quantities. We 
expect that improved data will allow us to determine the 
local value of the ratio r]/s. 

CONCLUSIONS AND OUTLOOK 

There are a number of issues that need to be addressed 
before an accurate value of 7] /s with fully controlled er- 
rors can be obtained. The most important of these is a 
better description of the transition from nearly perfect 
fluid dynamics in the center of the cloud to kinetic be- 
havior in the dilute corona. 

On the the theoretical side we need to develop tools 
that will allow us to perform calculations of transport 
properties in the strongly coupled regime. Some steps in 
this direction have been taken. Taylor and Randeiia de- 
rived sum rules for the viscosity spectral function [23]. 
Enss, Haussmann and Zwerger developed a resummed 
diagrammatic scheme that respects the sum rule and re- 
produces the kinetic limit [31]. There have also been 
some attempts at extending the AdS/CFT correspon- 
dence to non-relativistic conformally invariant theories, 
see [32, 33]. 



Finally, we would like to understand how nearly per- 
fect fluidity arises in different physical systems. In the 
case of the quark gluon plasma the question is whether 
nearly perfect fluidity is caused by strong interactions 
between well-defined quark and gluon quasi-particles, 
or whether the quasi-particle picture breaks down com- 
pletely and the low energy description involves non-local 
degrees of freedom, as in the AdS/CFT correspondence. 
In the case of the dilute Fermi gas we would like to un- 
derstand whether momentum transport is governed by 
fermionic quasi-particles, by collective modes (like the 
roton in liquid Helium), or whether there is no quasi- 
particle description at all. This question is difficult to ad- 
dress experimentally, but it can be studied numerically, 
by computing the viscosity spectral function. These cal- 
culations are still in their infancy, but we expect signifi- 
cant progress in the near future. 
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